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ABSTRACT 

Using information on the cosmic X-ray background and the cumulative light of active 
galactic nuclei (AGN) at infrared wavelengths, the estimated local mass density of 
galactic massive black holes (MBHs) and published AGN composite spectra in the 
optical, UV and X-ray, we compute the characteristic angular-integrated, broad-band 
spectral energy distribution of the average quasar in the universe. We demonstrate 
that the radiation from such sources can photoionize and Compton heat the plasma 
surrounding them up to an equilibrium Compton temperature (lb) of 2 x 10 K. 
It is shown that circumnuclear obscuration cannot significantly affect the net gas 
Compton heating and cooling rates, so that the above Tq value is approximately 
characteristic of both obscured and unobscured quasars. This temperature is above 
typical gas temperatures in elliptical galaxies and just above the virial temperatures of 
giant ellipticals. The general results of this work can be used for accurate calculations 
of the feedback effect of MBHs on both their immediate environs and the more distant 
interstellar medium of their host galaxies. 

Key words: galaxies: active - quasars: general. 



1 INTRODUCTION 

Massive black holes (MBHs) at cosmological distances came 
to our attention more than three decades ago due to the 
enormous outflow from them of high energy radiation. This 
observed fact, combined with the inverse square law, im- 
plies that the im mediate environs of these MBHs undergo 
dramatic heating jLevich fc Sunvaevlll97lt) when they are 
in the luminous "on" state - radiating at rates approaching 
Eddington limit for their masses. Yet most work on active 
galactic nuclei (AGN) or the environment of MBHs has ig- 
nored this " feedback" effect, even though it is easy to demon- 
strat e (e.g. ICowie. Ostriker fc Starklll97ll iPark fc Ostrikerl 
l200ll) that it should dramatically alter both the immediate 
environment of the MBH, from which accretion is occur- 
ring, and the more distant interstellar medium of the galaxy 
within which the MBH resides. The multiplicity of MHBs in 
the cores of galaxies found in groups and clusters may even 
sensibly alter the entropy floor of the ambient gas in these 
assemblages. 

But an accurate calculation of all of these effects could 
not be made until we had accumulated and averaged the 
spectral output of representative samples of AGN - and that 
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data has been missing; the lack can now be remedied due 
to recent observational advances. Furthemore, some means 
needed to be found to perform the averages, due to strong 
variability of the sources with time, with viewing angle and 
from object to object. 

It is the purpose of this paper to address these issues 
quantitatively: to compute the characteristic spectrum and 
all the important "Compton temperature", 2c, of the av- 
erage source and to estimate the consequences of exposing 
gas of cosmic chemical composition to the radiation from 
such sources. Specific applications of the general results of 
this paper will be presented in future publications. Jumping 
ahead to our conclusions, we constrain Tc to a narrow range 
around 2 x 10 7 K. This estimate is based on 1) measure- 
ments of the cumulative AGN light at various wavelengths, 
2) the measured local mass density of MBHs and 3) pub- 
lished composite quasar spectra. Interestingly, our estimate 
of the characteristic quasar Compton temperature is just 
above the virial temperatures of giant elliptical galaxies and 
somewhat below typical intracluster medium temperatures. 

Before turning to our detailed calculation of quasar 
spectral output, it may be useful to present some addi- 
tional background evidence and motivation. Comparison of 
the mass density of MBHs residing at the center of nearby 
galaxies with the total radiation flux from all AGN sug- 
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gests that MBHs have grown by radiatively efficient ac- 
cretion when their host galaxies were quasars. Overall ra- 



diation efficiencies [e. 



AJ?t/( AmBHC )1 of ~ 0.1-0.3 



are estimated ([Fabian fc I wasawd|l999fr [Salucci et al. Il999t 
lElvis. Risaliti fc Zamoranl2002l:lYu fc ^e^nainel2003) . It is 



important to note that this high mean efficiency is weighted 
by accretion rate; it is possible that during most of the life- 
time MBHs accrete much less efficiently, but most of the 
mass must be accreted in (perhaps brief) high efficiency 
episodes. 

Luminous quasars were much more abundant in the 
past (at redshifts ~ 2) than at the current epoch and one 
may reasonably be puzzled as to why the most massive lo- 
cal ellipticals (such as M87), which contain the most massive 
black holes, are so quiescent. Elliptical galaxies do contain 
significant amounts of hot (~ 10 7 K) gas capable of accreting 
onto central MBHs due to its fast radiative cooling compared 
to th e Hubble time. It has been suggested jBinnev fc Taborl 
Il995l) that feedback from the MBH may regulate accretion 
from the ambient gas, leading to an oscillation-type behavior 
such that short periods of strong nuclear activity are inter- 
changed with much longer quiescence periods during which 
the majority of observed galaxies are caught. This picture of 
brief intervals of high efficiency accretion separated by long 
periods of low level and low efficiency accretion is attractive 
in helping us understa nd the statistics of quasars and MBHs 
<Yu fc Tremainell2002l) . 

Continuing this theme, ICiotti fc OstrikeJ <1997t 1200 if) 
considered a scenario in which the gas of a central cool- 
ing flow is heated, during a quasar phase, above the galactic 
virial temperature by hard X-ray and gamma radiation from 
the MBH, which leads to a degassing of the central regions of 
the galaxy and switching off of the nucleus. A new episode of 
nuclear activity can begin after a large amount of cooled gas 
has again accumulated in the central regions of the galaxy. 
For this model to work, the characteristic quasar Compton 
temperature must be higher than the temperature of the 
cooling flow gas, otherwise the gas will be Compton cooled 
by the low-frequency radiation from the nucleus instead of 
being Compton heated. Regardin g the Compton tempera- 
ture as essentially a free parameter, ICiotti fc Ostriked (l200ll) 
presented solutions for a number of Tc values ranging from 
5 x 10 7 to 10 9 K, and the results are qualitatively indepen- 
dent of Tc so long ELS Cq = kTc I m v > vl , the stellar velocity 
dispersion. 

We note that the problem of gas preheating by X-rays 
emergent from MBHs finds its analogy in other astrophysical 
situation s: near stellar-mass b lack holes located in globular 
clusters dOstriker et al.| [l976) as well as in X-ray binaries 
with stell ar wind a ccretion onto a black hole or a neutron 
star JSunvaeJl97ll . 



2 BASIC ASSUMPTIONS AND 
CONSIDERATIONS 

2.1 Obscured vs. unobscured AGN 

In the standard AGN unification picture jAntonuccilll993l) . 
the active galactic nucleus is surrounded by an axisymmetric 
region (hereafter called the torus) filled with cold and dense 
material that intercepts and redistributes in wavelength and 



direction a substantial fraction of the primary radiation. The 
source will be classified as an unobscured (type 1) or ob- 
scured (type 2) AGN if its nucleus is observed directly - 
at a small angle to the axis of symmetry, or through the 
torus, respectively. One can also imagine a different situa- 
tion in which there are two physically different populations 
of sources: naked AGN (type 1) and AGN enshrouded in 
dense material (type 2). What plays the predominant role 
in nature, orientation or the presence/absence of a dense en- 
velope, is still a matter of debate, but in either case type 1 
and type 2 AGN exhibit distinctly different spectral energy 
distributions when observed far from the sources. 

It turns out that although circumnuclear absorption 
definitely plays a crucial role in shaping the observed spec- 
tral energy distribution of quasar emission, its effect on the 
characteristic Compton temperature is expected to be small. 
To demonstrate this let us cast the usual definition of the 
Compton temperature (e.g. iLevich fc Sunvaevlll97ll) - the 
plasma temperature at which net energy exchange by Comp- 
ton scattering between photons and electrons vanishes - in 
the following form: 



kT i rr VEF - dE +sr ^)EF Ed E 

4 J 10 ^F Ed E + J- kcV b { E)F E dE ■ 



(1) 



Here E is the photon energy, Fe(E) is the radiation spec- 
tral flux density, and the factors a(E) and b(E) represent 
Klein-Nishina corrections that become of importance in the 
hard X-ray regime (the corresponding explicit expressions 
are given in Appendix). Note that equation Q is valid in 
the limit kTc <§C m c 2 = 511 keV; it becomes inaccurate by 
more than 5% at kTc > 10 keV (in general the expression 
for Compton energy exc hange is nonlinear with re spect to 
gas temperature, see e.g. lSazonov fc Sunvaevl|200lf) . 

As will be detailed in and 21 f° r the characteris- 
tic spectral output of a type 1 quasar the Compton heating 
rate is completely dominated by the high-energy integral 
ZiokeV a {E)EFE dE while most of the Compton cooling is 
due to the component eV Fe dE. The latter integral rep- 
resents energy flux integrated over the 'blue bump' and the 
infrared band. As a result, we may expect the estimate 



kT c « - 
4 



lf^v a ( E ) EF z dE 



r 



Fe dE 



(2) 



to be accurate to a few per cent for type 1 AGN. 

If we now consider lines of sight passing through an ob- 
scuring torus, radiation emitted at E ;> 10 keV will be little 
affected by photoabsorption unless the torus is substantially 
Compton thick: Nn > a few 10 24 cm -2 . Further, most of 
the radiation emitted at E <J 10 keV will be absorbed and 
reemitted in the infrared, approximately retaining the total 
radiation flux due to energy conservation. Thus, neither the 
numerator nor the denominator of the RHS of equation (|5J 
is affected by obscuration to first order and we may expect 
that equation 10 with Fe representing as before the charac- 
teristic type 1 spectrum will also be a good approximation 
for the characteristic Compton temperature of type 2 AGN. 

The above line of argument is directly applicable to the 
scenario with two populations of sources, the unobscured 
and the obscured. On the other hand, orientation-based uni- 
fication schemes assume that an obscuring torus is present in 
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all AGN but it covers a solid angle of less than 4ir. In the case 
of such geometry, observers in unobscured directions will re- 
ceive both direct emission from the nucleus and a similar flux 
of reprocessed infrared radiation from the torus. As a result, 
the characteristic Compton temperatures of type 1 and type 
2 AGN could be somewhat lower and higher, respectively, 
than the angular averaged Tc given by equation (5) . 

We conclude that the characteristic Compton heating 
and cooling rates (per particle) should be the same within 
a factor of ~ 2 in all directions at a given distance from 
the AGN if the spectral distribution of primary emission is 
isotropic. This implies that knowledge of the angular aver- 
aged spectral output of a typical quasar would make it pos- 
sible to calculate with good precision the effect of Compton 
heating/cooling of gas around quasars. We note, however, 
that other radiative mechanisms such as photoionization 
heating and line cooling do depend strongly on the degree 
of circumnuclear absorption (see H5.2H . 

On the other hand, it is reasonable to expect that the 
two main spectral components of the intrinsic AGN radia- 
tion, the blue bump and the hard X-ray component, have 
somewhat different angular distributions. This will of course 
influence the ratio of the heating and cooling rates in equa- 
tion @. Since one of the main goals of this work is to obtain 
an unbiased estimate of the characteristic Compton temper- 
ature of the average quasar, we will give more weight in our 
computations to data on cumulative AGN light, such as the 
cosmic X-ray background and the total infrared flux from 
all AGN, than to another valuable source of information - 
composite spectra of type 1 quasars. The latter probably 
represent a snapshot of quasar emission in a certain cone of 
angles or/and from a certain population of objects. 

2.2 Cumulative AGN light 

As mentioned above, our derivation of the spectral output 
of the average quasar (in will be primarily based on 
measurements of cumulative AGN light from the sky. It will 
be based on the standard procedure described below. 

Let 6e{E,z) (erg s _1 cm -3 keV -1 ) be the angular in- 
tegrated spectral emissivity of AGN in a unit comoving vol- 
ume of the universe at redshift z. Then AGN located in 
the redshift interval [z,z + dz] will contribute to the locally 
measured surface brightness of the sky the amount 



dI B {E) 



47T 



dz 



(erg s 1 cm 2 sr 1 keV 1 ), (3) 



where t is the cosmic time and E' — E(l + z) is the energy 
a detected photon had when it was emitted. 

Integrating equation J5J over redshift gives the intensity 
of cumulative AGN light: 



Ie(E) = 



£e((1 + z)E,z)dz 



^HoJ (i + z)[n M (i + zy i + n A ] 



1/2' 



(4) 



We adopt the conco rdance model cosmology 
JOstriker fc Steinhardtlll995D : M = 0.3 and fi A = 0.7. 
Assuming that there is no spectral evolution of AGN with 
redshift, we may further write 



e B (E,z) = (F E (E))e(z). 



(5) 



Here, the function e(z) describes the evolution with redshift 
of the AGN comoving volume emissivity. 



Equations Q and establish the relation between 
the measurable spectral energy distribution of the cumula- 
tive AGN light, Ie{E), and the angular integrated spectral 
output of the average quasar, (Fe{E)), that we wish to de- 
termine and which characterizes the global energy release 
via accretion onto MBHs. 

The volume emissivity o f type 1 quasars reaches a 
maxi mum at z max ~ 2 (e.g. ISchmidt. Schneider fc Gunnl 
1995). Therefore, if the cumulative spectrum Ie{E) were 
known, one could roughly estimate (Fe{E)) simply by shift- 
ing Ie(E) to higher energies: E — * (l+z max )E. A more accu- 
rate treatment must be based on equations @ an d © and 
requires the knowledge of the evolution function e(z). We 
adopted the following evolution law in our computations: 



e(z) = 



(1 + *) 3 , 
e(zo) exp (zo 



z < z 
z > zo 



2.5. 



We have also considered the alternative law 



e(z) 



(1 + *) 
e(«o), 



z < z 

Z > Zo , 



(6) 



(7) 



We shall see (in iJ3.1ll that both types of evolution lead to 
practically the same results that are almost insensitive to 
the value of the critical redshift if zo > 2. 

The above parametrization is motivated by the follow- 
ing data. The evolution of type 1 quasa rs detected in the 
2dF and Sloan D igital Sky optical surveys (iBovle et all2000l : 
iFan et ai1l200ll) approximately matches the cutoff pattern, 
equation El, with z = 2.5. ROS A T soft X-ray (0.5-2 keV, 
iMivaii. Hasinger fc Schmiddl2000ft an d BeppoSAX hard X- 
rav f2-10 keV. lLa Franca et al.ll2002l) serveys suggest that 
type 1 quasars evolve approximately following equation J7J 
with zo — 1.5. The blazars observed at gamma-rays (above 
20 MeV) with CGRO/EGRET evolve as e(z) oc (1 + zf up 
to at least z ~ 2, with no d ata available for higher redshift 
l|Chiang fc Mukheriedll99gl) . 

The most recent Chandra deep X-ray (0.5-2 keV) and 
hard X-ray (2-8 keV) observations have confirmed a decline 
in the number density of quasars at z > 3 but addition- 
ally revealed a second peak in sou rce counts at z = 0.5-1 
jHasingerll200l ICowie et al]l2003ft . This low-redshift peak 
is, however, dominated by Seyfert galaxies, both unob- 
scured and obscured, with X-ray luminosities in the range 
Lx = 10 42 -10 44 erg s _1 , as compared to the previously 
known z ~ 2 maximum in the number density of quasars 
with Lx > 10 44 erg s _1 . For this reason, the low-redshift 
peak is less prominent in terms of the contribution to the 
total AGN X-ray light - in fact about a third of the cos- 
mic X-ra y background is no w believed to come from AGN 
at z < 1 jBarger et al]l2002>) . This could move the redshift 
of a typical AGN contributing to the cumulative AGN light 
from (z) ~ 1.5 (see H3.1|l to (z) m 1 but unlikely to lower red- 
shift. Since the relevant redshift correction factor is 1 + (z), 
the inferred average quasar spectrum could change by some 
25%. Considering this correction unimportant in view of the 
remaining uncertainties on the observational side, we ignore 
the z w 0.7 peak altogether in our calculations. 
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3 THE SPECTRAL OUTPUT OF THE 
AVERAGE QUASAR 

In this section we construct the average quasar spectrum 
from its pieces corresponding to different energy bands. 



3.1 High-energy emission 

We first consider emission above 2 keV. This spectral region 
is the main contributor to the Compton heating by quasars. 
It is now well established that the bulk of the cosmic X- 
ray and gamma-ray background (hereafter CXGB), the ob- 
served spectrum of which is shown in Fig. is composed of 
contributions from AGN of va rying degree of obscuration, 
luminosity and radio loudness j]VIadau 1 Ghisellini fc Fabianl 
11994 IComastri et alJll995t iGilli. Salvati fc Hasingerll200lli . 
This conclusion is based on 1) X-ray source counts, 
with some 80% of the E < 5 ke V background having 
been resolved into discrete sources jHasinger et alJ 1200 it 
Rosa ti et al ] 120021) . 2) data on the hard X-ray spectra of 
Seyfert galaxies, 3) the measur ed distribution of absorbing 
column densities in Seyferts iRisaliti. Maiolino fc Salvatil 



1999), 4) the redshift evolution of quasars and 5) for the 
gamma-ray part tentatively also on the observed similar- 
ity of the slopes of the CXGB spectrum and the spectra of 
EGRET blazars (ISreekumar et al]ll99St) . 

This allows us to identify, with some reservations as 
noted below (in &U.21 and J3.1.31 . the CXGB with the 
cumulative AGN spectrum Ie(E) and estimate the high- 
energy part of the average quasar spectrum (Fe) using the 
equations of H2.2I 

The CXGB spectrum was reliably measured in the 3- 
400 keV and 2 MeV-100 GeV bands by the HEAO-1 and 
CGRO observatories, and a usef ul formula fitting the mea- 
sured intensities was proposed bv lGruber et alJ l[l999). The 
spectrum peaks (when plotted in units EIe, see Fig. at 
~ 30 keV, is characterized by an approximately constant 
slope (Ie oc E~ a ) a w 0.3 at < 20 keV, by a varying slope 
q ~ 1.5 at energies above the peak up to ~ 10 MeV and is 
approximately a power law with a w 1.1 at higher energies. 

Extrapolating the Gruber et al. formula for the CXGB 
spectrum down to 1 keV, we find Ie (1 keV) = 8 keV s _1 
cm -2 sr _1 keV -1 , which is 20% below but is marginally 
consistent with the newest estimate based on all- sky obser- 
vations with RXTE/PCA (iRevnivtsev et all2003T) and simi- 
larly below the value derived from a joint analy sis of ROSAT 
and A SCA observations of small celestial fields jMivaii et alJ 
1998), the latter estimate being affected by cosmic variance. 
We thus apply the Gruber et al. fitting formula to describing 
the CXGB spectrum at all energies above 1 keV, except for 
the poorly explored 0.4-2 MeV region which we omit from 
our analysis, and adopt a 20% uncertainty in the CXGB 
normalization. 

We next adopt the template quasar spectrum 



{Fe) 

A 
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E l 
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B(l + kE -~<) E- 







2 keV <E<E 
E>E ={/3- a)Ei 
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/(I + kEg 



(8) 




The spectrum E{Fe) is shown in Fig. it peaks at ap- 



Figure 1. Measured CXGB spectrum (heavy solid line with a 
gap) and the adopted, rest-frame high-energy spectrum of the 
average quasar, equation JSJ (solid line). The redshiftcd light from 
a population of such quasars distributed in redshift according to 
equation JHJ has a spectral distribution that fits well the CXGB 
spectrum. For comparison shown are two other source templates, 
which generate similarly good fits to the CXGB spectrum when 
assuming different evolution scenarios: equation JBJ with zo = 3 
(dashed line) and equation J3 with zq = 2.5 (dash-dotted line). 



proximately 60 keV. The convolution of this average spec- 
trum with the evolution law given by equation © provides 
an excellent fit to the CXGB spectrum. We note that the 
characteristic slopes of the template J§) in different spectral 
regions are nearly the same as the corresponding original val- 
ues for the CXGB spectrum. This reflects the fact that only 
templates of the type ® can lead to the observed CXGB 
spectrum; essentially only the position of the peak is found 
from the fit. The high-energy template © is normalized so 
that (7<£;}(2keV) = 1. Our choice for the lower boundary 
of the spectrum (2 keV) is motivated by the expectation 
that obscured sources contribute negligibly to the average 
quasar spectrum at <; 2 keV and also by the fact that the 
high-energy component and the UV-soft X-ray component 
that will be discussed in H3.2l ioin near 2 keV in type 1 quasar 
spectra (Laor et al. 1997). 

Repeating the same fitting procedure for zq — 3 or as- 
suming the evolution given by equation @ with zo = 2.5, we 
obtain average spectra that deviate from the previous one 
by less than 10% (see Fig.0. This demonstrates a very weak 
dependence of the result on the assumed evolution scenario 
at z > 2.5, which reflects the fact that a typical quasar con- 
tributing to our average spectrum is located at (z) = 1.3-1.8 
(depending on E). Emission from substantially higher red- 
shifts contributes little to the background. 
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3.1.1 The average spectrum as a sum of the spectra of 
obscured and unobscured sources 

One can see that the average quasar spectrum just derived 
differs notably below a few tens keV from observed X-ray 
spectra of both unobscured and obscured AGN. This is, how- 
ever, consistent with this spectrum being the result of av- 
eraging over sources with a range of line-of-sight photoab- 
sorption columns. To illustrate this point and also for the 
purposes of our further radiative energy exchange computa- 
tions, let us introduce two additional spectra of which one 
will be representative of the spectra of type 1 quasars and 
the other will characterize their obscured counterparts. We 
impose the requirement that the appropriately weighted sum 
of these two spectra should match well the globally averaged 
one. We exclude from the current consideration blazars, as- 
suming that they do not contribute significantly to the av- 
erage spectra below ~ 500 keV, as will be justified in H3.1.2I 

Most of the published information on AGN X-ray spec- 
tra pertains to nearby (z <C 0.1), relatively low-luminosity 
(Lx = 10 42 -10 44 erg s _1 , the unabsorbed luminosity in 
the 2-10 keV ban d) Seyfert 1 and Seyfert 2 galaxies. 
IPerol a et alJ l)2002F) have recently presented a sample of 
broad-band (0.1-200 keV) spectra obtained with the Bep- 
poSAX satellite for nine Seyfert Is with Lx ranging from 
5 x 10 42 to 10 44 erg s . The measured spectral con- 
tinua above a few keV are well fitted by a model consist- 
ing of a power law with high-energy exponential cutoff, 
E~ a exp(— E/Ef), and a Compton reflection component. 
Typical values for the power-law index are a w 0.8 and 
the cutoff energies Ef span from about 70 keV to more than 
300 keV, typically E s ~ 200 keV. The relative amplitude of 
the reflection component attributed to a weakly ionized ac- 
cretion disk is R = 0.5-1. The above values of Ef are some- 
what lower than the previous estimates (Ef > 250 keV) 
based on EXOSAT GIN GA, HEAO-1 and CGRO/OSSE 
data JGondek et al.lll99rj) . We note that a rollover at ~ 50- 
100 keV was first clearly detecte d in the spectrum of 
NGC 4151 with GRAiVAT/SIGMA Jjourdain et alJll992lb 

A separate well-studied class of AGN is radio galax- 
ies, which can be considered the radio-loud counter- 
parts of Seyfert galaxies. Their X-ray spectra are sim- 
ilar to those of Seyferts, i.e. a m 0.8, Ef > 
100 keV, except that the reflection component seems 
to be weaker and is actually undetected in most cases 
jWc»znial£_el ; _alj_Jr^ fc Mushotzkvl 

l2000l: iGrandi. Urrv fc Maraschill2002lb Note that the stud- 
ied radio galaxies are on average more luminous (Lx ~ a 
few 10 44 erg s -1 ) than the Seyferts discussed above. 

On the other hand, the available information on lumi- 
nous (L x £ 10 44 erg s" 1 ) AGN located at z > 0.1, which 
usually are classified as quasars, is rather scarce and gen- 
erally limited to E < 20 keV. However, the AGN X-ray 
luminosity function is characterized by an increasing typ- 
ical luminosity L* with redshift, so that L+(z < 1) ~ a 
few xlO 43 erg s" 1 while L*(z > 1) ~ a few xlO 44 e rg s _1 
jMivaii. Hasinger fc Schmid J2000I: ICowie et al.ll2003i) . As a 
consequence, a typical AGN contributing to the CXGB is a 
quasar at z ~1— 2 with a luminosity Lx J> 10 44 erg s _1 . 

All of the available information on quasars (exclud- 
ing blazars) is consistent with their X-ray spectral con- 
tinua being very similar to those of local Seyfert and ra- 
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Figure 2. Heavy solid line: Adopted high-energy spectrum of the 
average quasar, equation JSJ, with the gap at 1-5 MeV approx- 
imately corresponding to the 0.4—2 MeV poorly explored region 
in the CXGB spectrum (see Fig.0. Dashed line: Adopted spec- 
trum of the average unobscured quasar, equation J5J. Dotted line: 
Adopted spectrum of the average obscured quasar, equations 1101 
and 1111 , Note that above ~ 50 keV the obscured and unob- 
scured spectra coincide because photoabsorption becomes neg- 
ligible. Solid line: Weighted sum, equation 1121 . of the average 
obscured and unobscured spectra. Dash-dotted line: Estimated 
upper limit, equation 1131 . on the contribution of blazar emission 
to the average quasar spectrum. 

dio galaxies. In particular, the 1-20 keV spectra of both 
radio-quiet and radio-loud type 1 quasars are characterized 
by a ss 0.8 jLawson fc Turnerlll997j;lReeves fc TurnedEoOol : 
iHasenkopf. Smabrriiil^^E^cleou^l2002Tl ~In addition, high- 
energy spectral cutoffs with Ef ~ 100 keV have been re- 
ported for at least two nearby (z ~ 0.1) quasars, PG 1416- 
129 a nd MR 2251-178 iStaubert fc MaisadJl99dlOrr et ail 
l200ll) . It appears that the Compton reflection component is 
generally weak at l east in the most lumino us quasars with 
L x ~ 10 46 erg s _1 jReeves fc Turnerll200Clb 

Summarizing the above facts, we adopt the following 
spectrum (plotted in Fig. |2J for the average unobscured 
quasar: 

/El \Ci H 771-O.8 -B/200 keV , n \ 

(F B )(typel) oc E e 1 (9) 

The chosen value for the cutoff energy not only falls 
comfortably in the range of Ef values found for nearby AGN 
and quasars, but also leads to the right location of the peak 
in our globally averaged quasar spectrum. We also note that 
the addition of a R ~ 0.5 reflection component to equa- 
tion © would have a fairly similar effect on the spectrum 
as slightly modifying the main parameters — am 0.7 and 
Ef ~ 150 keV. Since such changes are within the current 
observational uncertainties and of the same order as source- 
to-source variations, we ignore the reflection component. We 
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also ignore the usually observed fluorescence K Q iron line, 
since its equivalent width is small (~ 100 eV, the above 
references) . 

We next introduce another spectrum, which we at- 
tribute to the average obscured quasar: 



(F B )(typel) 

x I f(N B )exp[-o-(E)N H ]dN H , (10) 



<F B )(type2) 



where o-(E) represents the photoelectric absorp- 

tion cross section for sola r chemical composition 

^Morrison fc McCammo We adopt a lognor- 

mal distribution of absorption columns centered at 
N H = 10 24 cm" 2 : 



(11) 



/(log JV H ) = 4= ex P [-OogJVH - 24) 2 ] . 

Our choice of the central value for /(logA%) is in first 
place dictated by the shape of the globally averaged quasar 
spectrum at E < 50 keV, which is much flatter than the 
unabsorbed spectrum In addition, the median value 
of the inferred Ah distribu tion for the local p opulation of 
Seyfert 2s is ~ 10 24 cm -2 jRisaliti et allll999h . The width 
of the Ah distribution, equation I1H . was taken just suffi- 
ciently large to produce a relatively smooth spectrum. For 
simplicity, we neglect here the Compton scattering of nu- 
clear radiation in the obscuring torus, which can signifi- 
cantly modify the emergent spectru m when log A?h <; 24.5 
jMatt. Pompilio fc La Francalll999l) . 

The absorbed spectrum given by equations HOH 
and < I U is shown in Fig. We should mention that 
X-ray spectra exhibiting a low-energy exponential cut- 
off have been recently measured for several obscured 
quasars (iwasawa. Fabian fc Et tor1 l200ll: Is tern et al.l l2002t 
iNorman et al.ll2002l) . 

Finally, we demonstrate in Fig. [5] that our globally av- 
eraged quasar spectrum is well fit below ~ 300 keV by a 
weighted sum of the adopted average spectra of type 1 and 
type 2 quasars: 



(F E ) w 0.25<F B )(typel) + 0.75<F B )(type2). 



(12) 



This corresponds to the ratio 3 : 1 of obscured to unob- 
scured sources, similar to that (4 : 1) estimated for the local 
population of Seyfert galaxies llMaiolino fc Riekdll995|) . 

3.1.2 Contribution from blazars 

Although there is little doubt that in the 2-300 keV band our 
adopted average quasar spectrum represents mostly normal, 
nonblazar quasar emission, the situation is more uncertain in 
the gamma-ray region where blazars and non-AGN sources 
may contribute to the average spectrum inferred from the 
CXGB. 

According to a popular view fe.g. lZdziarsklll996fl . the 
CXGB above a few MeV is to a large part due to blazars, 
an AGN class which we excluded from our consideration in 
H3.1.1I The spectra of blazars are believed to be dominated 
by Doppler boosted radiation fro m a relativistic jet point - 
ing close to our line of sight (e.g. lUrrv fc Padovariilll995l) . 
as compared to the quasi-isotropic emission from normal 
quasars. Regardless of whether powerful jet emission is pro- 
duced in most quasars but is usually beamed away from us, 



or it occurs only in a small fraction of them, obviously it is 
important to estimate the contribution of blazar emission to 
the average quasar spectrum. As we shall demonstrate in 
using a constraint obtained below, blazar beamed emission 
is at best a minor contributor to the global AGN energy 
output as well as to the global radiative heating by AGN. 

In contrast to normal quasars, it proves impossible to 
produce any meaningful average spectrum for blazars, since 
blazar spectra vary dramatically from object to object and 
also with time: they usually consist of two broad peaks, 
whose positions and relative amplit udes vary by orders of 
magnitude fe.g. iFossati et al]ll99Sl) . Nonetheless, we can 
still derive an interesting upper limit on the expected con- 
tribution from blazars to the average quasar spectrum using 
the following facts: 

(i) The 30 MeV-10 GeV spectra of some 50 blazars de- 
tected by CGRO/EGRET are consistent with being simple 
power laws, with the a verage energy index a m 1.15 ± 0.04 
( Muk heriee et alJll997f) , which is equal within the errors to 
the slope (a « 1.10 ± 0. 03) of the CGRB spect rum in the 
30 MeV-100 GeV range JSreekumar et al]ll99SD . 

(ii) An estimated luminosity function of gamma-ray 
blazars and its evolution, combined with the spectral in- 
formation above, suggest that blazars contribute at least 
a signi ficant fraction of the CXGB in the EGRET ener gy 
range iChiang fc Mukheriedri998USreekumar et allll998ft . 

(iii) The spectra of gamma-ray blazars peak at Me V en- 
ergies (Schonfcldcr 1994; Mc Naron-Brown et alJ ll995). 

(iv) The spectral slopes of gamma-ray blazars mea- 
sured between 2 keV and a few x 10 2 keV vary consid- 
erably fro m object to object but tend to group around 
a sa 0.7 ^McNaron-Brown et alJ Il995l : iKubo et al] Il99gl : 
iTavecchio et alJl2002l) ~ 

(v) The contribution of all types of radio-loud quasars 
to the CXGB backgrou nd in the ~ 1 keV range is <. 10% 
idella Ceca et all I1994T) . and the contribution of blazars, 
thoug h very uncertain, is likely sever al times lower, ~ 1% 
llComastri. Di Girolamo fc Settil ll996 ) . 

The following spectrum (see Fig.|2Jl therefore provides a 
plausible upper limit on the contribution of blazar emission 
to the average quasar spectrum: 

,„ uw s / 0.082£-°- 7 e- B/1 ° 4kcV , 2keV<£'<4MeV 
(JfiMblaz = < . , ' — 

x /v J \ 1.52S- 11 , £>4MeV. 

The spectrum 113H was obtained by fixing the lower-energy 
and high-energy slopes at a = 0.7 and 1.1, respectively, 
and also assuming that the blazar contribution to the glob- 
ally averaged spectrum is 5% at 2 keV and ~ 100% at 
> 100 MeV. These conditions uniquely determine the char- 
acteristic break energy. 



3.1.3 Contamination by other types of sources 

The largest uncertainty is associated with the 300 keV- 
10 MeV part of the average spectrum (see Fig. First, 
we essentially do not know the intensities at 1-5 MeV (see 
Fig-HJ since the CXGB has not been reliably detected in the 
corresponding 400 keV-2 MeV band (here we assume that 
a typical source contributing to the CXGB is at z ~ 1.5). 
Moreover, other types of astrophysical sources, not neces- 
sarily AGN, may provide the major contribution in this 
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energy range. The most pro mising candidate proposed so 
far is Type la supernovae llThe . Leising fc Clavtonlll993t 
iRuiz-Larmende. Casse fc Vangioni^iamrl2n0 J) . Finally, it 
cannot be ruled out that a subclass of blazars exhibiting 
a peculiar spectral bump at ~ 1 MeV may provide a signif- 
icant contribution ijBlom et alJll99fll : IComastri et alJll99r^ . 
The last hypothesis remains highly speculative at present, 
since only a few such "MeV blazars" have been observed. 

Fortunately, although the uncertainty associated with 
the average quasar spectrum above 300 keV is large, it can- 
not affect significantly our estimates of the global heating 
and cooling rates by AGN (sec [J1J. 



3.2 Near-infrared to soft X-ray emission 

We next consider the 1 eV-2 keV band (which corresponds 
to wavelengths A from 1.2 /im to 6 A). Since there are prac- 
tically no data on cumulative AGN light for this spectral 
range, our approach below will be different than in the pre- 
ceding section. 

In the standard AGN unification scenario, an obscuring 
torus of dust and gas with Ah ~ 10 22 -10 24 cm~ 2 will be 
transparent to X-rays above ~ 10 keV and also to infrared 
and low-frequency radiation at A J> 10 /im. All the near- 
infrared to soft X-ray radiation emitted by the nucleus along 
obscured lines of sight will be absorbed and reemitted at IR 
and submillimeter wavelengths. We may thus expect that 
the spectrum of the average quasar at 1 eV-2 keV will be 
completely dominated by emission from unobscured sources. 
We can therefore take advantage of the ample existing in- 
formation on type 1 quasar spectra, which we may condense 
in the simple spectral form 



(F E ) = 1.20 



159£T 



1 eV < E < 10 eV 
10 eV < E < 2 keV. 



(14) 



Here 10 eV approximately corresponds to Xh ya = 1216 A 
and the adopted normalization is such that the combined 
spectrum given by equations © and 1141 is continuous at 
2 keV. 

In deriving equation 1141 . we took into account miscel- 
laneous recently published data on quasars, i ncluding the 
Hubb le Space Telescope catalog of UV spectra jTelfer et alJ 
2002), composite optical-UV spectra fro m the Sloan Digital 
Sky S urvey (SDSS, Vanden Berk et al. IVanden Berk et alJ 
l200ll) . 2-2 keV spectra from the ROSAT Bright Quasar 
Survey llLaor et al.l Il997ft . and statistics on the character- 
istic 2500 A-2 keV spectral slope , aox (lYuan et alJll998t 
IVignali. Brandt fc Schneiderll20p3|). We also used the atlas 
of quasar energy distributions o flElvis et all <ll994f) . We note 
that the combined contribution of resolved emission lines, 
including Lya, to the total luminosity of the blue bum p 
dominated by the continuum is small jTelfer et al]l2002ft . 
but nevertheless an attempt has been made to take it into 
account in the definition 1141 . 

The aox index is particularly relevant for our study, 
since it together with the spectrum of the hard X-ray com- 
ponent determines the ratio of the amplitudes of the hard 
X-ray and blue bumps in the quasar spectru m. Equation 
1141 yields aox = 1.4. For comparison, the lElvis et alJ 
(1994) mean radio-quiet and radio- loud quasar spectra are 
characterized by aox = 1.4 and aox = 1.3, respectively. 



From analysis of optical and ROSAT, Chandra and XMM- 
Newton X-ray data for optically sel e cted SDSS quasars 
(with redshifts up to 6~). IVignali et al.l <|2003T) have inferred 
a dependence of aox on quasar rest-frame UV luminosity: 
10 30 erg s" 1 Hz" 1 and w 1.7 for 



1.5 for L 



2500 A 

"'Hz 



aox 

L ' ~ lO 32 erg s _i Hz _i . We note that these monochro- 
2500 A 

matic luminosities correspond to L2-iokeV ranging from 
~ 10 44 to ~ 10 46 erg s _1 and, as was noted in H3.ll the lower 
values from this range are characteristic of the quasars pro- 
ducing the bulk of the global AGN energy output, suggesting 
that aox ~ 1.5 might be a typical value globally. 

We note that the template given by equation 1141 pos- 
sesses a number of key properties of observed quasar spec- 
tra, including a turnover near ALya llTelfer et all2002T) and a 
gradual flattening in the 0.2-2 keV band. The adopted tem- 
plate, primarily inferred from recent UV observations, is red- 
der than the 'big blue bump' predicted by most theoretical 
accretion disk models whose spectra ext end to the extreme 
ultraviolet (see iKoratkar fc Blaezl 1 19991 for a review) . We 
should also note that the spectra of radio-quiet and radio- 
loud quasars are quite similar up to ~ 100 eV felvis et alJ 
I1994T) . with « 90% of all quasars being radio quiet and this 
fraction being see mingly independent of redshift and b olo- 
metric l uminosi ty ijSteni^^lj boOOl : llvezic et alJl20o3 : see 
however ICirasuoioet^ni2003Tr 

Given the remaining uncertainty in the aox value char- 
acterizing the globally averaged quasar spectral output, we 
estimate that the ratio of the optical-UV to high-energy in- 
tegrated radiation fluxes that follows from equations J8J and 
1141 for the average quasar spectrum is uncertain by a factor 
of ~ 2. 

We note that recent data seem to indicate that the scat- 
ter in the near-IR to soft X-ray spectrum from quasar to 
quasar is fairly small. In deed, according to t he Hubble Space 
Telescope observations llTelfer et alJ 1200*21) the blue bump 
peaks invariably at A ~ 1200 A (i.e. not far from 2500 A, the 
wavelength appearing in the definition of aox) for quasars 
whose monochromatic luminosities vL v at 1100 A range 
from a few 10 45 to a few 10 47 erg s _1 (which corresponds to 
£2-10 kcv ranging from ~ 10 44 to ~ 10 46 erg s" 1 ). On the 
other hand, the rms scatter in apx for a given luminosity is 
of the order of 0.1 l|Vignali et alJ l2003'l. which corresponds 
to a factor of 2 scatter in the ratio i^L„(2 keV)/^L„(1200 A), 
which is similar to our estimated uncertainty in this ratio 
for the average spectrum. 



3.3 Medium-infrared to submillimeter emission 

Finally, we consider the spectral range A > 1 ^m, a major 
contributor to the net Compton cooling. As for the high- 
energy region, both obscured and unobscured sources are 
expected to significantly contribute to the integrated AGN 
light at IR and submm wavelengths. Although there is only 
fragmentary information on Ie(E) in this case, we shall use 
it below to derive some constraints on the average quasar 
spectrum at A > 1 /im. There remains, after all available 
observations are utilized, a substantial uncertainty in the 
medium-IR luminosity of a typical quasar. However, as will 
be noted and demonstrated below, we have an integral con- 
straint that provides an additional strong bound. The total 
luminosity of all quasars is limited by the observed, low- 
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redshift density of MBHs. This limits the integrated low- 
frequency emission. 

Our approach will be similar to several previous studies 
of the contribution of AGN to the cosmic IR and submm 
background (hereafter CIB. [ Almaini. Lawrence fc Bovld 
ll999tlRisaIiti. Elvis fc Gillill2002h . We take into account the 
following data and limits related to the cumulative IR quasar 
light: 

(i) The CIB spectrum (after subtraction of the cosmic 
microwave background) at 125-2000 fim mea sured by the 
COBE/FIRAS experiment jFixsen et al.lll998l) . 

(ii) Luminous AGN, which produce the bulk of the cos- 
mic X-ray background, contribute less than 10% to the CIB 
background at 850 /im, as inferred from the cross-correlation 
of SCUBA submm sources with Chandra deep X-ray sur veys 
JSevergnini et alJl200d : iFabian. Smail fc Iwasawall2000l) . 

(in) The CIB at 15 urn: EIe (15 urn, total) sa 3 nW m~ 2 
sr -1 flFranceschini et al 1 120011: lElbaz. Cesarskv fc Chanlail 
120021) . and « 17% of this intensity is due to AGN that 
make up ~ 85% of the 2-10 keV background, as inferred 
from the cross-correlation of XMM-Newton and Chan- 
dra deep X-ray surveys with ISOCAM infrared surveys 
jFadda. Flores fc Hasingerl2002l) . One can thus estimate the 
cumulative quasar mid-IR light as EIe(15 |im) = 0.5 ± 
0.15 nW m" 2 sr" 1 jFadda et alJl2002l : lElbaz et al.ll2002l) . 

(iv) Upper limits on the CIB at A < 40 fim set 
by TeV cosmic opacit y measurements toward blazars 
iFranceschini et ai1l200ll) . 

The above data are displayed in Fig. [3] We also show 
the expected cumulative quasar light spectrum in the near- 
infrared (A < 1 /im), calculated from the blue bump tem- 
plate, equation 1141 . with the normalization fixed by the 
cosmic X-ray background. Evidently only loose constraints 
can be obtained on the far-IR emission of the average quasar 
from the presented data alone. Fortunately, the radiative en- 
ergy exchange between the quasar and the surrounding gas 
does not depend at all on the spectral shape at IR wave- 
lengths. The only important thing is the total luminosity 
emitted in this spectral band, and it is possible to obtain 
some constraints on this quantity from the observational 
data presented in Fig. [3] 

Indeed, it is widely accepted that the mechanism re- 
sponsible for the IR emission of radio-quiet quasars (and 
partly for that of radio- loud q uasars) is therm al radiation 
from radiatively heated dust l|Barvainisl Therefore, 
whichever is the spectral distribution of the IR emission 
from a quasar or from an ensemble of quasars, it cannot 
be narrower than that of modified black body emission from 
monotemperature dust: 

E 3+P 

FE(dust) oc 



cxp(S/fcT d ) 



(15) 



where T d is the dust temperature, /3 ~ 1-2 is the emissivity 
index, and we have assumed the emission to be optically thin 
[the effect of a significant optical depth can be mimicked by 
varying the values of j3 and Td in equation 115^ ]. 

We show in Fig.[3]several examples of spectra of the cu- 
mulative quasar light computed from equations (|1J|-(|SJ for 
the template spectrum l|15p taking different Tj values. Nor- 
malizing the resulting spectra at 15 /im to 0.5 nW m -2 sr" 1 , 
the current best estimate for EIe, we can derive an upper 
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Figure 3. Observational constraints on the spectrum of cumu- 
lative AGN IR and submm light: the data point with error bars 
at 15 fim and the upper limit at 850 /an. Also shown are the 
COBE/FIRAS spectrum of the total extragalactic background 
(CIB) at 125-2000 (im (thick solid line) and upper limits on the 
CIB at 3-40 )ira. See text for references to the data. All these data 
cobined allow a large variety of possible IR spectra of the aver- 
age quasar, including the modified black body emission spectrum, 
equation 1151 : the result of convolution of this model spectrum 
with the quasar evolution function, equations J1J— JH}, is shown 
for P = 2 and three values of the dust effective temperature (as 
indicated in K above the corresponding curves) that either max- 
imize or minimize the integral EIe- The additional integral con- 
straint given by equation 1221 leads to the adopted average quasar 
IR spectrum, equation 1231 (dashed line, after convolution with 
the quasar evolution function). Also shown (dotted line) is the 
adopted average quasar spectrum in the blue bump region, equa- 
tion 1141 . after convolution with the quasar evolution function, 
with the normalization fixed by the cosmic X-ray background. 



limit on the integrated IR emission from quasars as a func- 
tion of the characteristic dust temperature Td. In this way we 
find that J I E dE < 15 nW m~ 2 sr" 1 and < 8 nW nT 2 sr" 1 
for Ti = 100 and 1500 K, respectively. These limiting dust 
temperatures lead to maximal luminosities. The more gen- 
eral upper limit J Ie dE <J 20 nW m 2 sr 1 is applicable to 
any realistic spectral distribution, characterized by a range 
of T d . 

On the other hand, there is a more stringent, lower limit 
on the cumulative, frequency-integrated far-IR light from 
quasars: 



Im.= 



I E dE > 0.8 ± 0.3 nW m" sr~ 



(16) 



This minimum corresponds to dust emission with Td ~ 
400 K (see Fig. El , and the quoted error results from the 
uncertainty in the observed flux of cumulative AGN light at 
15 /im. 
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Let us now derive an upper limit on /ir, following the 
argument of lSoltar] lll982T) . The mass density of MBHs in the 
local universe is estimated as pbh = (3± 1) x 10 5 Mp, Mpc -3 
<Yu fc Tremaineil2002l: lAller fc Richstoneil20oil . taking h = 
0.72 ± 0.05 iSpereel et alJ l2003h . Assuming that all of 
this mass has been accumulated through accretion with a 
redshift-independent mass-to-radiation conversion efficiency 
e 7 gives an upper limit on the total, frequency-integrated ra- 
diation flux from all quasars: 

e 7 c 3 p B H / e(z)/(l + z)\dt/dz\dz 



I E dE < 



4tt 



(1.7 ±0.6) 



J e(z)\dt I dz\dz 
nW m~ 2 sr"\ 



(17) 



where we have again used equation Q to describe the quasar 
emissivity evolution (although the dependence of the result 
on e(z) is weak). 

Now, the total flux /total consists of contributions from 
the high-energy region (> 1 keV), the blue-bump region 
(1 eV < E < 1 keV) and the low- frequency region (< 1 eV). 
The first of these components can be found directly by in- 
tegrating the CXGB spectrum shown in Fig. 



iX-ray 



(18) 



where the quoted error has been propagated from the un- 
certainty in the CXGB normalization. 

A more conservative estimate can be obtained by in- 
tegrating the X-ray background over the 1-150 keV range, 
where it is certainly dominated by AGN emission: 



h 



-isokev = 0.21 ± 0.04 nW m~ 2 sr~\ 



(19) 



The corresponding integral over the blue bump can be 
found by substituting the template given by equations JSJ 
and 1141 1 into equations and using the normalization 

provided by the CXGB: 



/bine = 0.4±^ nW nr 2 sr-\ 



(20) 



where we have roughly estimated the observational uncer- 
tainty in the average optical to soft-X-ray spectrum of type 
1 quasars. 

We can now estimate the maximum possible integrated 
quasar infrared (A > 1 mm) light: 



IlB. — /total — /blue — II- 150 keV 

< (1.7 ± 0.6) -0.7 ±0.3 nW m~ 2 sr~\ 
~ v y 0.1 



(21) 



One can see that the upper limit on /ir, equation l)21|l . 
is consistent with the lower limit on /ir, equation 1161 . if 
e 7 is above 0.05, implying that MBHs have grown by ra- 
diatively efficient accretion via a standard disk. A similar 
conclusion was e arlier reached by oth er researchers (e.g. 
lElvis et al]l2002t lYu fc Tremaindl2002ft . If we additionally 
require that e 7 < 0.1, then we find from the limits given by 
equations 1161 and 1211 that 



/ir = (1.0 ±0.6) nW m~ 2 sr~ 



(22) 



It was implicit in the above treatment that the ra- 
diative output of quasars is purely the result of accretion 
onto MBHs. Recent observational data lend support to this 
view, seemingly excluding starburst activity as the domi- 
nant contributor to the bolometric luminosity of powerful 
AGN producing the bulk of the cosmic X-ray background. 



Indeed, high-quality spectra obtained wi th the ISO satel- 
lite for dozens of optically bright quasars dHaas et alJl200d : 
IPolletta et"aill200fil : lAndreani et alJl2003h . plotted in units 
EFe, exhibit a broad bump rising from 1-2 /im (followed 
by the blue bump on the short-wavelength side) and de- 
clining in most cases at A > 60 fim. This bump is inter- 
preted as multitemperature thermal emission from warm 
dust, with the maximum temperature (1000-1500 K) rep- 
resenting the evaporation temperature of dust grains. In 
addition, the far-IR spectra measured for s everal obscured 
quasars also tend to peak below 30 /im (e.g. lBarvainis et alJ 
ll995tlDeane fc Trentharrll200ll) . 

The fact that the bulk of the infrared luminosity of 
quasars is emitted at wavelengths shorter than 60 (jm im- 
plies that dust cooler than ~ 50 K is typically not ener- 
getic ally important. This, a ccording to the standard theory 
fe.g. lRowan-Robinsonlll995Tl . suggests that the observed in- 
frared emission is mostly the result of reprocessing by dust 
of the optical-UV radiation from the nucleus, rather than 
resulting from starburst activity. 

Overall, the available data are consistent with the dom- 
inant contribution to the infrared luminosity of a typical 
quasar being due to dust heated to 50-1500 K by the radia- 
tion released as a result of accretion onto the central MBH. 
For this reason we consider equation 1221 a robust estimate 
that can be used to normalize the average IR quasar spec- 
trum, for which we adopt the following template: 



(F E ) = 1.5xl0 4 



£ 5 T d " 7 



K exp(£/fcT d ) - 1 



dT d , E < 1 eV, (23) 



where E and Td are measured in keV and K, respectively. 

This spectrum, which bears a greater resemblance to 
observed quasar spectra than equation 1151 . represents op- 
tically thin modified (/3 = 2) black-body emission of dust 
characterized by a range of temperatures (100-1400 K) and 
constant emitted luminosity per unit logarithmic Td interval 
(so that the spectrum has a flat core when plotted in EFe 
units). The numerical coefficient in equation 1231 . which is 
defined relative to the X-ray component, equation JHJ, en- 
ables that f o lcV I E dE = 1.0 nW m" 2 sr" 1 . The adopted 
upper Td value enables continuity at 1 eV between the dis- 
tributions given by equations 1231 and 1141 . while the lower 
T d boundary and the spectral form have been chosen rather 
arbitrarily to allow the predicted cumulative quasar light to 
be consistent with the data point at 15 ^im (see Fig. 



3.4 The broad-band spectrum 

We are now finally in a position to put together the different 
pieces of the average quasar spectrum. The resulting broad- 
band spectrum is shown in Fig. [I] It has three broad maxima 
in terms of energy output: one in the medium infrared, the 
blue bump and another one at ~ 60 keV. It is necessary to 
note that the radio (A > 1 cm) component of the quasar 
spectrum is unimportant, because its contribution to the 
total energy output and consequently to the net Compton 
cooling rate is very small. This follows from the estimated 
surface brightness of the extr agalactic radio background - 
f I E dE ~ 10" 4 nW m" 2 sr" 1 <Bridldjl967HDwek fc Barker! 
2002) and also from spectroscopy of individual quasars - 
the radio luminosity of radio-quiet and radio-loud quasars 
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Figure 4. Solid line: Broad-band spectrum of the average quasar 
adopted in this paper. Dashed line: Adopted spectrum of the av- 
erage unobscured quasar, rescaled by a factor of 1.1 for visibility. 
Dotted line: Adopted spectrum of the average obscured quasar, 
rescaled by a factor of 0.9. The weighted sum of the adopted 
type 1 and type 2 spectra (not shown) matches well the average 
spectrum below 300 keV. 

is typically 3 and 5 orders of magnitude smaller, r espectively, 
than the luminosity of the IR or blue bumps (e.g. lElvis et al. 
Il994l) . We note, however, that in the vicinity of parsec-scale 
quasar jets simulated Compt on heating by their radio radi- 
ation can become important dLevich fc Sunvaevlll97ll) . 

Since we shall study in the following sections the ef- 
fect of AGN obscuration on gas heating and cooling, we 
show in Fig. 2] two additional spectra that may be consid- 
ered representative of type 1 and type 2 quasars. We have 
defined these templates in a somewhat arbitrary manner, 
which is unavoidable for several reasons, in particular the 
distribution of emitted energy between the hard X-ray com- 
ponent and the blue bump may be intrinsically anisotropic, 
i.e. depend on the line of sight even in the absence of cir- 
cunmuclear absorption; also the obscuration pattern may be 
more complicated than assumed, e.g. some quasars may ap- 
pear obscured in the UV but unabsorbed in the X-ray, etc. 
Specifically we required that 

• the ratio of obscured to unobscured quasars be 3 to 1, 

• in the X-ray region, the type 1 and type 2 spectra be 
given by equation @ and equation 1101 . respectively, 

• the type 2 spectrum have no emission in the blue bump 
region, 

• the infrared bump be the same for all three considered 
templates. 

The weighted sum of the adopted obscured and unobscured 
quasar spectra matches well the globally averaged quasar 
spectrum below ~ 300 keV. We recall that this upper bound- 



ary results directly from the cutoff energy (200 keV) that we 
adopted in equation © and that there remains a significant 
observational uncertainty in the position of this cutoff for 
Seyfert galaxies and quasars. 



4 COMPTON TEMPERATURE 

We can now determine the Compton temperature Tc char- 
acterizing the spectral output of the average quasar, which 
is the gas temperature at which Compton heating balances 
Compton cooling. An accurate result in the limit kTc >C 
m B c? — 511 keV can be obtained from equations given in 
Appendix by requiring that (dW/dt) + + (dW/dt)- = 0. 

Let us note here that if additionally E <C m e c 2 , 
the Compton cooling rate becomes proportional to the 
frequency-integrated flux J {FE)dE. We may thus expect 
that for our adopted average quasar spectrum most of the 
cooling will be provided by the infrared component, with ad- 
ditional smaller contributions coming from the optical-UV 
and high-energy bands. Klein-Nishina corrections act to fur- 
ther reduce the cooling by hard X-rays and gamma-rays. On 
the other hand, the Compton heating rate is, in the limit 
E <C m c c 2 , proportional to another integral, J E{FE)dE, 
which is completely dominated by the high-energy compo- 
nent of the quasar spectrum. However, also in this case 
Klein-Nishina corrections dramatically diminish the heat- 
ing by photons with E > 100 keV. 

Both the Compton heating and the Compton cooling in- 
tegrals are well constrained for the adopted average quasar 
spectrum. First, the spectral shape of the high-energy com- 
ponent is known very well at least up to 300 keV and its 
normalization is uncertain to within ~ 20% - see equation 
118H . Second, the uncertainty in the radiation flux integrated 
over the entire average spectrum is ~ 40%, as follows from 
equations 1161 . 1171 and 121H . We can thus determine the 
characteristic Compton temperature quite accurately: 

T c = (1.9 ±0.8) x 10 7 K. (24) 

We emphasize that the surprisingly small uncertainty 
in the final result has resulted from three facts: 1) good 
knowledge of both the spectrum and normalization of the 
CXGB, 2) the lower limit on the cumulative AGN light in 
the infrared band that follows from the cross-correlation of 
deep X-ray and infrared surveys and 3) the upper limit on 
the total flux from all quasars inferred from the local den- 
sity of MBHs. It is interesting that the detailed properties of 
the blue bump play a minor role in obtaining equation 124H . 
Indeed, the contribution of the blue bump to the Comp- 
ton heating rate is negiligible, while its contribution to the 
Compton cooling rate is of the order of 25%, which is less 
than the uncertainty in the Tc value for the average quasar. 

The value given by equation 12 H corresponds to the en- 
tire spectrum extending into the gamma-ray band. If we ex- 
clude from consideration the spectral segment above 300 keV 
whose AGN origin still needs to be proven or disproven, we 
find 

T c (< 300 keV) = 1.3 x 10 7 K. (25) 

Therefore, AGN emission above 300 keV contributes less 
than 30% to the total Compton heating rate (its contri- 
bution to the Compton cooling rate being negligible). The 
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contribution of emission above 500 keV is < 20% and this 
excludes blazar beamed emission as a major contributor to 
the global heating due to AGN. 

We can similarly calculate the Compton temperatures 
corresponding to the adopted type 1 and type 2 quasar spec- 
tra (see Fig. 0J: 

0.8 x 10 7 K, 



Tc(unobsc) 
Tc(obsc) = 2.0 x 10 7 K. 



(26) 
(27) 



Because we have taken substantial freedom in defining the 
type 1 and type 2 spectra, the above values are clearly more 
uncertain than the Compton temperature corresponding to 
the globally averaged spectrum. 

We should note that the above calculation of Compton 
energy exchange is valid in the case of a Maxwellian mo- 
mentum distribution of electrons, and thus assumes that the 
time scale for isotropization of electron momenta by mutual 
collisions is shorter than that for scattering of hard pho- 
tons by electrons. We also explicitly assumed in our calcula- 
tions that there is no significant bulk motion of the plasma 
caused by the gravity and radiation pressure of the central 
MBH. Both of these assumptions can be violated very close 
to quasar nuclei. 



5 EFFECTS OF EXPOSING GAS TO THE 
RADIATION FROM THE AVERAGE 
QUASAR 

Having built the spectrum and determined the characteristic 
Compton temperature of the average quasar, we can now es- 
timate the consequences of exposing gas of cosmic chemical 
composition to the radiation from such sources. 

5.1 Maximum distance of heating in the low 
density limit 

As was noted in ^ it is well established that MBHs sitting 
in the centers of local galaxies have experienced an epoch or 
multiple epochs of rapid growth accompanied by the appear- 
ance of a quasar. One may therefore ask: to which maximum 
radius could a MBH currently of mass Mbh heat during its 
life time gas from a low initial temperature to the quasar 
Compton temperature? Such formulation assumes that the 
gas is of sufficicnetly low density that it is fully photoionized 
so that only Compton heating and cooling are important. 

For gas of temperature T the heating rate per electron 
at a distance r from the quasar is 

dW _ L(t)a T 4fc(T c - T) 



4:Tvr 2 



m c c z 



(28) 



where L is the quasar luminosity, and we have assumed that 
the source is isotropic. Since the quasar is powered by ac- 
cretion, the luminosity can be related to the rate of growth 
of the MBH: 



r /+\ dM&K M 2 



(29) 



7 dt 

Hence, the total energy received during the growth of the 
MBH by an electron-proton pair located at r is 



ATI/ kTc OT 2 
AW = - rfc,MBHC , 



assuming that T < Tc. 

If we now require that each electron-proton pair receive 
an energy at least AW = 3kTc, we find the maximum dis- 
tance out to which this can be done: 

a T e 7 MBH\ 1/2 . ( ^ V' 2 f M } x ' 1 



rc 



(H^L) 1/2 = o.4kpc(^V 



(31) 



5.2 Heating/cooling of a partially ionized gas 

Up to now our attention has been focused on Compton heat- 
ing and cooling. However, gas exposed to intense quasar ra- 
diation may be sufficiently dense to remain only partially 
ionized. Photoionization heating as well as cooling through 
continuum and line emission will then be important. We 
consider this situation below. 

We assume that the gas is optically thin and is in ion- 
ization equilibrium. Note that the characteristic ionization 
and recombination times are typically much shorter than the 
time scale for Compton h eating/cooling. We adopt the sola r 
element abundances from lGrevesse. Noels fc Sauvall il996f) . 
unless stated otherwise. Under these assumptions, the ion- 
ization balance as well as the heating and cooling rates are 
fully determined by the instantaneous gas temperature, the 
radiat ion spectral distribution an d the ionization parameter 
jTarter. Tucker fc Salpeterlll969h . which we define here as 



i±£ =i.4x lO 9 ^- 



Mbh 1 cm 



1O S M n 



(32) 



9 9 f 



(30) 



where Lboi (erg s" 1 ) is the bolometric, angular-integrated 
luminosity of the central source, Z/Edd is the Eddington lu- 
minosity for a given MBH mass, n the hydrogen nucleus 
density and r the distance from the source. To proceed we 
assume that 3/4 of the whole sky as seen from the MBH is 
obscured while 1/4 is clear and that the emergent spectrum 
is described by either the adopted type 2 or type 1 spectrum 
(see Fig. |3}, dependent on from which part of the sky the 
source is observed. According to this picture, observers lo- 
cated at the same distance but in different directions from 
the source will receive equal radiation fluxes above ~ 20 keV 
and below ~ 1 eV, but those looking through the obscuring 
torus will not see any optical, UV or soft X-ray radiation. 

We performed comp utations using the latest version of 
XSTAR (iKallmanl 120021) . to which we have added a block 
responsible for calculating Compton heating and cooling 
from equations 1A11 and IA2t . Our code also accounts 
for an additional heating caused by Compton scattering 
of hard X-rays off electrons bound in hydrogen and he- 
lium atoms. To this end we use the simple approximation 
that the Compton heating rate per electron is the same 
for bound and free electrons and is described by equa- 
tion JAlft . This implies that the net Compton heating rate 
per unit volume does not depend on the fraction of free 
electrons in the plasma, as opposed to the Compton cooling 
rate, which is proportional to this fraction. This has been 
shown to be a good ap proximation for photons with ener- 
gies above a few keV feasko. Sunvaev fc Titarchukl Il974l : 
ounvacv fc Churazovl ll996). and such photons completely 
dominate the Compton heating in our case. 

Since the definition of £ adopted in XSTAR is different 
from ours, equation 1321 . we provide here the relations be- 
tween Lboi and the ionizing flux F 1QTL (between 13.6 eV and 
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Figure 5. (a) Dependence on the ionization parameter, 
equation 1321 . of the stationary temperature of plasma of 
solar chemical composition exposed to the different radia- 
tion spectral distributions considered in this paper: the glob- 
ally averaged quasar spectrum (solid line), type 1 spec- 
trum (dashed line), type 2 spectrum (dotted line), and addi- 
tionally (F B )(type2) + 0.01(T B )(typel) (dot-dashed line) and 
(F B )(type2)+0.05(F B )(typel) (short-dash-long-dashed line), (b) 
The effect of plasma chemical composition on the T B t a t(£) curve 
for the average quasar spectrum. Solid line: solar element abun- 
dances, dotted line: iron is absent, dashed line: H-He gas. 



13.6 keV) in obscured and unobscured directions as well as 
its average over the sky: 



4nr 2 F io 



(47vr 2 F ion ) 
Lhoi 

■ (unobsc) 



= 0.13, 



0.50. 



4nr 2 F io 



-(obsc) = 0.012. (33) 



Suppose now that the irradiated gas has had enough 
time to achieve thermal equilibrium. Fig. |2K shows the gas 
stationary temperature T sta t as a function of £ for each of 
the different spectral templates presented in Fig. 2] Fig. 
illustrates the effect of gas enrichment by metals and iron on 
the r at at (C) curve. We see that the gas reaches the Compton 
temperature given by equations I24I - 127H when £ > 10 . 
In this region the gas is fully photoionized and the ther- 



mal balance is dominated by Compton heating and cooling. 
At £^ 10 5 , the Ttat (5) 

curves are fairly similar for the 
globally averaged spectrum and for the type 1 spectrum, 
and are mainly determined by the balance between pho- 
toionization heating and various cooling mechanisms (see 
iKallman fc McCravll982l for a detailed discussion of the un- 
derlying physics) . 

Therefore, for given luminosity Lboi and gas radial dis- 
tribution n(r), r(£ = 10 5 ) defines the size of a Compton 
heating zone where gas can achieve a steady state with 
T = Tc. Consider as a specific example M87, a giant el- 
liptical galaxy in the Virgo cluster. M87 hosts a 3 x 1O 9 M0 
MBH and contains hot (T ~ 10 7 K) gas characterized by 
n ~ 0.1, 0. 05 and 0.02 cm~^ at r = 1, 4 and 10 kpc, re- 
spectively dMatsushita et aT1l2002l) . For these parameters, 
the outer boundary of the Compton heating zone would be 
at r(10 5 ) w 1 kpc if the M87 nucleus switched on at its 
Eddington luminosity. This size is somewhat smaller than 
rc = 2.2 kpc given by equation 131H for M87. 

The Tstat (0 curve is quite different for £ < 10 5 in the 
case of the type 2 spectrum, for almost complete lack of 
ionizing UV and soft X-ray photons. In this case, there 
is a narrow transition region at £ ~ 10 dividing a high- 
temperature and low-temperature regions characterized by 
Tstat = T c ~ 10 7 and ~ 10 4 K, respectively. Moreover, 
the solu tion in this transition zone is known to be un sta- 
ble (e.g. lBuff fc McCravHl974l : lKallman fc McCravl 1983) . It 
is interesting to note that in the low-temperature region, 
where the gas is practically neutral, Compton heating due 
to scattering on bound electrons contributes 8% to the to- 
tal heating rate dominated by photoionization of metallic 
atoms by hard X-rays. This, however, has negligible effect 
(less than 1%) on the resulting equilibrium temperature in 
the low-temperature region. The effect is even smaller for 
the average and type 1 quasar spectra, in which case the 
photoionization heating associated with the radiation above 
13.6 eV is several orders of magnitude stronger than Comp- 
ton heating on bound electrons in the low-ionization limit 

(£<io- 4 ). 

The results presented in Fig. [KJi thus imply that radia- 
tive heating can be quite different for £ <J 10 5 , i.e. outside 
the Compton heating zone, if plasma is exposed to the ra- 
diation reprocessed by the obscuring torus instead of being 
irradiated directly by the quasar nucleus. The direct conse- 
quence of this is the possibility of significantly anisotropic 
radiative feedback on the quasar environement in the sce- 
nario with partial obscuration of the nucleus: strong pho- 
toionization heating (of an initially cold material) can take 
place within the cones of direct nuclear emission but not in 
obscured directions. 

We, however, consider the above situation unlikely, be- 
cause it resulted from the complete absence of UV and soft 
X-ray radiation in our adopted type 2 spectrum. In reality 
the heating anisotropy will be significantly reduced if a small 
fraction of the ionizing radiation outgoing from the nucleus 
within the open cones is scattered into obscured lines of 
sight. This is demonstrated by computations we performed 
for input spectra given by (i 7 B)(type2) + 0.01(T£;)(typel) 
and (FE)(type2) +0.05(Tfi)(typel). This corresponds to the 
scattering gas having a Thomson optical depth of tt ~ 0.04 
and ~ 0.2, respectively, assuming that this gas is ionized (so 
that electron scattering of UV radiation is possible) and fills 
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the central funnel of the torus and extends some distance 
above its top, covering at least / = 1/4 of the sky as seen 
from the nucleus. As can be seen from Fig. QJ,, the addi- 
tion of such a small amount of scattered ionizing radiation 
changes the situation dramatically, significantly diminishing 
the departure of the resulting T s tat(£) curve from those cor- 
responding to the globally averaged and type 1 spectra. In 
such a case, irradiation by the quasar will probably have 
rather similar effects on gas located along unobscured and 
obscured lines of sight not only inside the zone governed by 
Compton heating and cooling but also outside it. 

Different kinds of observation do reveal significant quan- 
tities of warm (T ~ 10 4 5 -10 6 K) ionized gas located at dis- 
tances from possibly a fraction of a pc up to hundreds of 
pc fr om the nuclei of Seyfert galaxies (e.g. iKrolik fc Kriisl 
2001). In Seyfert 2 galaxies such as NGC 1068, this plasma 
reveals itself via polarized scatt ered emission in the opti- 
cal (e.g. iMiller fc Goodrichlll990t) and via scattered contin- 
uum, recombination emission and resonance scattering i n 
the X-ray band (e.g. iTurner et al"]ll997l : lOgle et alJkoOSl) . 
The same photoionized gas produces absorption features in 
the UV and soft X-ray spectra of Seyfert 1 galaxies an d 
some quasars (e.g. iGeoree et al1ll99St iMathur et a!lll994) . 
Futher, ~ 10% of optically selected quasars show broad ab- 
sorption lines (BALs), arising in material apparently flow- 
ing outward from the nucleus with velocities ~ 10 4 km s , 
and it appears that mos t or possibly all q uasars contain 
such BAL outflows (e.g. iGreen et all 1200 ID . Finally, pow- 
erful radio galaxies at high redshifts show strong UV line 
and continuum emission extended along the axis of the ra- 
dio source, which is interpreted as scattered rad iation from 
an obscured quasar nucleus (e.g. lTran et alll998h . Taken to- 
gether, all these diverse observations suggest that a warm or 
hot gas characterized by the product Trf < a few per cent - 
of the same order as the scattering fractions assumed in our 
simulations above - may be ubiquitously present in AGN. 
We also note that electron scattering optical depths tt ~ 
several 10 -2 are predicted for the ~ 10 pc cooling flows 
feedin g quasars in elliptical galaxies (e.g. ICiotti fc Ostrikeil 
1200 ll) . 

We finally note that the T sta t(C) curves shown in Fig.0 
descend below 10 4 K when £ becomes less than ~ 1. This im- 
plies that the gas density must exceed 10 7 cm _3 (10 pc/r) 2 
at distance r from a 1O 8 M0 black hole radiating at Edding- 
ton limit in order to provide conditions for the existence of 
a dusty molecular torus such as postulated in the AGN uni- 
fication model. Such dense material cannot fill more than a 
small fraction of the obscuration region, otherwise its large 
Thomson optical depth tt > 2 x 10 (10 pc/r) will prevent 
even hard X-rays from escaping through the torus and make 
the source Compton thick. We note that XSTAR becomes 
unsuitable at T < 3000 K. 



6 CONCLUSIONS 

The main results obtained in this work are as follows. 

We combined information on the cumulative AGN light 
in the IR and X-ray bands, the estimate of the local mass 
density of MBHs and composite optical to soft X-ray quasar 
spectra to construct in a robust way the angular-integrated 
radiation spectrum of the average quasar in the universe 



( 35p . This spectrum characterizes the global energy release 
via accretion onto MBHs, and is the result of implicit sum- 
ming over unobscured and obscured sources. 

We calculated ([0Jl the Compton heating and cooling 
rates for gas exposed to radiation with the adopted average 
quasar spectral distribution. The Compton heating results 
from downscattering in energy of the hard X-rays, while the 
cooling is due to inverse Compton scattering of primarily 
the IR photons. The Compton temperature, representing the 
equilibrium between Compton heating and cooling, is well 
constrained near 2 x 10 7 K, with an estimated uncertainty 
in this value of ~ 50%. 

We presented simple arguments (Q and supported 
them with accurate calculations (S0J that circumnuclear 
obscuration cannot significantly affect the Compton heat- 
ing and cooling rates. As a result, the Compton tempera- 
tures characteristic of obscured and unobscured directions 
(or sources) are within a factor of 2 from the average Tc 
value quoted above. 

The almost invariant shape of the (unabsorbed) hard X- 
ray spectrum combined with the recently published data of 
observations of tens and hundreds of optically bright quasars 
that cover the blue bump segment of the spectrum up to its 
connection with the X-ray component suggest that the rms 
source-to-source scatter in the ratio of the amplitudes of the 
hard X-ray and blue bumps is fairly small, ~ 2, at least 
for type 1 quasars. Of the same order is the possible trend 
going from relatively low-luminosity (Lx ~ 10 44 erg s _1 ) 
to high-luminosity (Lx ~ 10 46 erg s _1 ) quasars ( H3.2I . We 
thus believe that the average spectrum and Compton tem- 
perature found in this work should also represent well the 
spectral output of individual quasars, at least of their major- 
ity. We note that the IR emission from quasars results from 
reprocessing of the primary UV emission, so that both com- 
ponents always have comparable powers, and therefore the 
variability of the shape of the IR spectral component from 
source to source has practically no effect on this conclusion. 

Blazar-type beamed emission is energetically unimpor- 
tant globally, contributing on average less than 20% to the 
net Compton heating rate f H3. 1.211 . Nevertheless, current ob- 
servations do not rule out the possibility that this additional 
radiation component can play an important role in gas heat- 
ing for particular objects. This question needs further study. 

We showed f H5.ll that during its lifetime a MBH of 
mass Mbh can heat to the Compton temperature low- 
density gas within a radius of ~ 0.5 kpc(MBH/lO 8 M0) 1//2 , 
which would normally constitute a significant fraction of the 
effective radius of the spheroid in which the MBH resides 
and a negligible fraction of the core radius of a cluster of 
galaxies. 

We performed computations ( H5.21 of the radiative 
heating of partially photoionized plasma of cosmic chem- 
ical composition exposed to the radiation from the aver- 
age quasar, taking into account photoionization heating and 
plasma cooling through line and continuum emission in addi- 
tion to Compton heating and cooling. Although the derived 
Tatat(C) curves are quite different for the adopted obscured 
and unobscured spectra, we demonstrated that scattering of 
a relatively small fraction (~ a few %) of the primary nuclear 
radiation by dense gas in the vicinity of the active nucleus 
will wipe away most of this difference. Observations indi- 
cate that such scattering does occur in most AGN. There- 
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fore, circumnuclear obscuration probably can be ignored to 
first order in considering the radiation feedback of MBHs on 
their environment. 

The general results of this work provide us with the 
basis for follow-up work aimed at accurately calculating the 
feedback effect of MBHs on their gaseous environment. 



APPENDIX A: 

In the limit kT <C m c c 2 — 511 keV, the characteristic rates 
per electron of plasma heating and cooli ng due to spon- 
taneous Compton scattering ar e given by dShestakov et all 
ll988l:ISazonov & SunvaelkOOll) 



/dW\ 
V dt ). 



(a;-3)(x + l)ln(2s + 1) 
dx 



€E{X) L8*3 
-10x 4 + 51a; 3 + 93a; 2 + 51x + 9 



4a; 2 (2a; + 1) 



and 
fdW\ 
V dt J . 



aTin c c 
-a-rckT 

€ E (x) 



5 



3(3x 2 - 4a; - 13) 



dx 



(Al) 



ln(2a; + 1) 



-216x b + 476a; 5 + 2066a; 4 + 2429a; 3 
8x 2 (2a; + l) 5 

1353a; 2 + 363a; + 39 

8a; 2 (2a; + l) 5 
4a T ckT 

47x 



dx 



+ 0(x 2 ) 



dx, 



(A2) 



where x = E/m c c 2 and £e(x) is the radiation spectral en- 
ergy density. It follows that a(E) = 1 — 21E/5m c c 2 + ... and 
b(E) = 1 - 47E/8m e c 2 + ... in equations Q and ©. 
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